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Abstract: In this paper we discuss the force exerted by the field of an 
optical cavity on a polarizable dipole. We show that the modification of 
the cavity modes due to interaction with the dipole significantly alters the 
properties of the force. In particular, all components of the force are found 
to be non-conservative, and cannot, therefore, be derived from a potential 
energy. We also suggest a simple generalization of the standard formulas for 
the optical force on the dipole, which reproduces the results of calculations 
based on the Maxwell stress tensor. 
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1. Introduction 

Starting with the pioneering works of Ashkin p|2| , where trapping and manipulation of objects 
by optical forces was first demonstrated, there has been an explosion of interest in optical 
forcing of particles. Optical tweezers have since been developed into a standard tool used in 
many applications, while on the fundamental side optical cooling of atoms has opened unique 
opportunities for exploring various quantum mechanical many-body phenomena. Recently a 
great deal of interest has been devoted to the possibility of optical cooling of macroscopic 
objects such as mechanical nanoresonators |(3j|6), membranes or particles oscillating in 
optical traps |T(F-T3^|. While traditional optical trapping experiments involve freely propagating 
laser beams, optical microcavities, which confine light in a small volume, have emerged as a 
candidate for the source of optical forces 1 16 17). The spatial confinement of light within the 



small mode volume of the cavity results in an increase of the strength of the force. Additionally, 
a material object interacts with the cavity mode, resulting in a shifted modal frequency and 
altered field distribution. As a result, the dynamics of mechanical degrees of freedom become 
coupled with those of the field giving rise to so-called dynamical backaction 1 18] responsible 
for effects such as cavity cooling or heating of mechanical degrees of freedom ||5, |12„19J . 

Understanding cavity optomechanical phenomena depends on correct representation of the 
optical force exerted by the cavity modes. In the case of free-propagating optical fields (i.e. 
laser beams), the force on a subwavelength object (dipole) is naturally separated into gradient 



and scattering components |20|. The gradient component is analogous to the force on a static 
dipole, which tends to draw a particle into regions of greater field intensity. It can be presented 
as the gradient of the electromagnetic energy of the polarized particle and is therefore conserva- 
tive. The scattering component results from radiation pressure and is expressed in terms of the 
momentum flux impinging on the particle per unit time. This force is non-conservative because 



it results from the process of irreversible exchange of momentum and energy between the par- 
ticle and optical field. Due to its conceptual simplicity and apparent universality, this paradigm 
has become firmly engrained in the current literature and has been accepted as a framework 
for calculating optical forces also due to cavity modes. The effects of spatial confinement are 
taken into account by using cavity modes to represent electromagnetic field while allowing the 
resonant frequencies of the cavity to depend on mechanical degrees of freedom. 

However, in our recent paper |21| we pointed out that the validity of the conserva- 
tive/gradient, non-conservative/scattering paradigm explicitly depends on the assumption that 
the particle itself does not change sources of the incident field. While this condition is fulfilled 
for free-propagating electromagnetic fields, it is violated for fields confined within a cavity. 
The position-dependence of the cavity resonance frequencies is just one manifestation of this 
phenomenon. In order to elucidate all consequences of the particle-induced modification of the 
cavity field, we considered the interaction between whispering-gallery-modes of a spherical mi- 
croresonator and a small dielectric particle. In this system forces can be calculated by a rigorous 
analytical approach based on the Maxwell stress tensor. These calculations show that the force 
cannot be described within the standard paradigm |21]. In particular, no vector component of 
the optical force is conservative, i.e. derivable as a gradient of potential energy. Furthermore, 
the force tangential to the surface of the cavity, which is responsible for the "carousel" effect 
observed in Ref. |22|, was found to have a contribution proportional to the real part of the 
particle's polarizability. This contradicts the assertion of Ref. p2{ that the tangential optical 
force, which is in the direction of the momentum flux of the unmodified by the particle cavity 
mode, is of strictly "scattering" origin. In Ref. |21 1 these results were obtained using formal 
ab initio approach, which while providing accurate expressions for the force, does not allow 
for the simple physical interpretation of the results and for finding connections between them 
and heuristic approaches used in works of other authors. The objective of this paper is to show 
how the traditional gradient-scattering force approach can be generalized to derive the results 
of Ref. f21 1 without relying on complicated Maxwell stress tensor calculations. It can be con- 
jectured that the "pseudo-gradient" formalism of this paper can be applied even in situations in 
which rigorous treatment is not possible. In addition, the results presented here allow elucidat- 
ing limitations of previous heuristic approaches to optical forces. 

To achieve this objective in the most efficient way, we organize this paper in the following 
manner. We begin by reviewing the derivations of the gradient and scattering forces paying 
particular attention to the assumptions involved. Based on this discussion we propose a psuedo- 
gradient procedure as a way to extend standard gradient/scattering approach to optical cavities. 
We use this idea to determine the force on a particle due to a spherical whispering gallery mode 
resonator and compare it with the results of stress tensor based calculations. 

2. Optical force on a small polarizable particle 

2.1. Gradient force: thermodynamic derivation 

We begin by recalling the thermodynamic approach to deriving the electrostatic force on a small 
dielectric particle. This approach emphasizes one of the key requirements for the validity of the 
gradient paradigm. Following the classical textbook by Landau & Lifshitz p3| the electrostatic 
component of the free energy of the polarizable particle is 



where D = CqE + P and the integral runs over the volume of the particle. E here refers to the 
total field in and around the particle, and P is its polarization. This total energy contains the 
energy, Ug^, of the external field, £, which would have existed in the absence of the particle. 




(1) 



Subtracting Ugxt, the energy required to polarize the particle is 



Upoi = u,ot -^J eo\(B\^dV 
Using this expression one can derive the change in Up^i due to a small variation of 2 as 

5up^i^- fp-5€dV. (2) 



The force is determined by substituting in 5upoi the change in the field, 5€, due to an infinites- 
imal translation of the particle, dvp, which is given by: 5£ = {drp ■ V)C;. Equation |2] then takes 
the following form: 

5up„i^-5rp- [(P-V)(BdV. (3) 



Taking into account that this change in energy can be related to the work of the force F acting 
on the particle: 5upoi = — F • 5rp, one obtains expression for the force as 



F = J{P-V)edV. (4) 

If £ is approximately constant over the dimensions of the particle (dipole approximation), it 
can be taken out of the integral, giving 

F = (p-V)£, (5) 

where the dipole moment p = /PdV. Assuming linear polarizability of the particle, p = ai£, 
one derives the final gradient expression for the force: 

F = ^aoV|e|2 (6) 

It should be noted, however, that this derivation depends critically on the assumption that a 
displacement of the particle does not affect the distribution of the sources of the external field 
Without this assumption one would not be able to equate the change of € due to small 
variations of the sources (Eq. (j2])) with its change due to particle's displacement making the 
subsequent equations invalid. 

2.2. Gradient force: direct derivation 

An alternative derivation, which is also commonly encountered in the textbooks, is based on a 
model of an electric dipole as a system of equal and opposite charges -izq, separated by some 
small, ultimately infinitesimal distance d. The same approach can be used to describe forces 
on an induced dipole characterized by polarizability a. The dipole is assumed to be placed 
in some external field E, so that p = aE. We make no demands on E other than that p be 
defined self consistently with it. In particular, E may be dependent on the dipole itself For 
example, if the external field E is due to charges on a conductor, the presence of the dipole 
will alter the charge distribution and thus E. We make this explicit by writing E = E(r,rp), 
where r and are respectively a field point and the position vector of the particle. The total 
force is derived by considering the Coulomb forces at each charge comprising the dipole: F = 
g'[E(rp-|-d/2,rp)-E(rp-d/2,rp)] (see Fig.[T|i. By taking the limit |d| ^0, keeping |p| = 
^|d| constant, the force derived in this case is: 



F=[(p.Vr)E],^r, 



(7) 



where refers to a gradient with respect to field coordinates r. While this result looks similar 
to Eq. (|5]l, there is an important difference between them. The force given by Eq. Q implies 
that before taking the spatial derivative of the field with respect to field coordinates r, one sets 
coordinate of the particle rp to coincide with r. On the other hand, Eq. (|7j) requires that this 
procedure is reversed. Physically, it reflects the fact that the electric force on a dipole results 
from the spatial variation of the electric field across it. The two equations, Eq. (jSj) and Eq. (j?]), 
produce identical results only if the field exerting the force does not depend upon particle's 
position. 

Equation (j7]i can also be transformed into a "gradient" form 

F=iaVr|E(r,rp)|2|r=r,. (8) 

which, however, differs from Eq. (j6]l. Unlike the latter, Eq. ([8]l involves taking the gradient 
of the function of two variables, and, therefore, the expression a|E(r,rp)p/2 cannot be inter- 
preted as a potential energy unless E is independent of Vp. Correspondingly, the force calculated 
according to Eq. (|6]l does not have to be conservative. 

To simplify terminology and notations in the subsequent consideration we will call the oper- 
ation presented in Eq. m a "pseudo-gradient" and will use notation V to represent it. 




Fig. 1 . Set up for evaluating the force on a dipole modeled as a system of equal and opposite 
charges. The distance between the charges will be taken to zero. 

2.3. Total force on a dipole 

\n this subsection we generalize the previous results to the case of a dipole interacting with a 
time-dependent harmonic electromagnetic field. The oscillating dipole moment creates a cur- 
rent density d^/dt, giving rise to a magnetic force. The total Lorentz force on the particle for 
incident fields independent of the dipole's position 03 = — // a)V x 2; have the standard form 
F = (p • V) € + lip/iif X 58 . The time averaged expression for this force can be rewritten as |24) : 

(F) = ^i^e[a]V|(£|2 + ]^Jm[a] {(0^e[€* x «B] + J^m[{€* • V)£]) . (9) 



where a now is a complex valued (with radiative corrections included) polarizability of 
the dipole. Based upon analysis of the previous sub-section we conjecture that the expres- 
sion for the force in the case of the field dependent on the particle's position {E,B} = 
{E(r,rp),B(r,rp)} can be obtained from Eq. (Ml by substituting 2;,*8, V with E,B, V respec- 
tively. For a standard dipole particle with radius Rp, refractive index n^, and polarizability 

2 

a = Aneo{aQ + -ik^al), 

where 

«0=^p(np-l)/(np+2), 
Equation ([9]l for the force can be re-written as 

(F) =-V(M) + f7c(g)+(7c:^^m[(E*-V)E], (10) 

2(0 

where (m) is the average polarization energy of the dipole (m) = — |^e[p-E], a is its scattering 
cross section a = ,f m[a\k / Eq, and (g) is the average momentum density of the field (g) — 
j^e[eoE X B], and c is the speed of light in vacuum. 



In the subsequent sections of the paper we will apply Eq. ( 10 1 to the case of a small dielectric 
particle interacting with a whispering gallery mode (WGM) of a spherical resonator Comparing 
the obtained expressions for the force with results of Ref. |21 1, where the optical force in this 
system was calculated using Maxwell stress tensor approach, we will be able to substantiate 
validity of Eq. ( [T0| and shed additional Ught on physical properties of the optical forces due to 
cavity-confined electromagnetic field. 



3. Optical force of a WGM resonator 

3.1. WGMs of a single spherical resonator 

An optical whispering gallery mode is a long living excitation of a spherical resonator, which 
can be thought of as a ray of light propagating along the equator of the sphere and trapped in 
it due to total internal reflection. In a spherical coordinate system with polar axis perpendicu- 
lar to the plane of the propagation of the mode {XYZ system in Fig. |2]), its field is described 
by a single vector spherical harmonic (VSH) characterized by polar, azimuthal, and radial in- 
dexes /,m,i respectively and polarization, TE or TM. For concreteness we focus only on TE 
polarized modes with radial index s= \ defined as E = Eo^h^/\kr)Xi ,„{9 , 0) (time harmonic 
factor e^""' is assumed and suppressed). Here £0 is a normalization factor, ^ is the scatte- 
ring amplitude describing response of the resonator to an incident radiation with frequency (0, 
/ij'^ (kr) is the spherical Hankel function of the first kind, and k is the magnitude of the wavevec- 
toi k — (o/c. In the close vicinity of a chosen resonant frequency the scattering amplitude can 
be approximated as 

^ = -/r{"V(«-«r+'r|''^) (11) 

with coj^^ and yj^^ being respectively the resonant frequency and decay rate of the mode. The 
angular portion of the VSH is defined as X/ „, = LF/ „,(0, 0)/y^/(7+Ty, where L is the di- 
mensionless angular momentum operator L = — /r x V, and Yi ,„{9,<p) is the scalar spherical 
harmonic. We focus here on so called fundamental modes with m = I. WGMs with long life- 
time are characterized by Z ^ 1 and in expressions that follow we neglect terms of order l/l 
compared to those or order unity. 

For calculation of optical forces due to this WGM it is also convenient to find an expression 
for its field in a coordinate system with polar axis connecting the center of the sphere and 



Fig. 2. Coordinate systems used for evaluation of optical forces together with schematic 
presentation of the resonator, WGM, and the dipole. The axis always connects the center of 
the resonator and the point of observation. When using this coordinate system to calculate 
the forces, this axis passes through the center of the particle. 



the point of observation {X'Y'Z' system in Fig. |2]) Such an expression can be obtained using 
rotational properties of VSH p5] expressed as 



m'=—l 



(12) 



where D|^)^^(a,)3,7) is the Wigner D function and a,j3,7 are the Euler angles specifying the 
rotation from the unprimed to primed coordinate system. The D-functions are defined by 



D^'] {a.B.y)^e 

w'm ^ ' r ' i / 



-i{m' a+mf) 



d" (J3), 



where the function dl^\ (/3 ) for m ~ I can be written as 



(2/)! 



{l + m)\{l-m)\ 



B B 

cos — sin — 
2 2 



cot(^ 



(13) 



In order to find an expression for the WGM in the primed system one has to apply inverse 



transformation 7, — j3,— a) ~ [D)^j'^^,(a,j3,7)]* to the VSH defined in the unprimed 

coordinate system. Rotation by the angle 7 is equivalent to shifting the 0' coordinate, so we 
only consider transformations whith 7=0. Applying this transformation to a WGM with orbital 
number I — Lwe find its representation in the rotated system in the following form: 



E' — Eoh^i^^Y^a^iXi, 



(14) 



where 



= ifD|^[(0,-j3,-a) = 



yo + ' 



(15) 



Here we introduced dimensionless variable yo = (o) — (»|"')/rj|'' representing the relative de- 
tuning of the external frequency from the resonance of the WGM with respect to its width. At 
any point on the polar Z' axis (9' — 0), all ,„ vanish except for Xl.±i- Assuming L ^ 1, and 
dispensing with the prime on E', the field can be written explicitly as: 



E: 



-.Eah^l\kr 




(16) 



where E,± — {iO ±(j))/ y/2 and 9 and are the spherical coordinate unit vectors referred to the 
global, unprimed system. The magnetic field H = B//io = — //o/ZqV x E is given by: 



H-£o 




(1)/ 



kr kr 



{kr)]' 



(17) 



where the prime denotes differentiation with respect to argument kr. We can express the field at 



any point in space in the form of Eq. ( 16 1 and Eq. ( 17 1 by changing the Euler angles appearing in 
a,„. It is important to note that angles a and j3 correspond to the respective angular coordinates 
(j) and 9 of the Z' axis as viewed from the unprimed coordinate system. 

3.2. Calculation of the force neglecting particle-resonator coupling 

In order to elucidate the effects of the particle-induced modification of WGMs on the properties 
of the optical force, we first compute this force with this modification neglected. In this case, 
assuming that a particle with coordinates Tp ~ {rp,9p,(j)p) as defined in the XYZ coordinate 
system, lies on the Z' axis of the X'Y'Z' system, we can substitute the results of the previous 



sub-section, Eq. (16i and (17i, into Eq. (lOi for the force, while replacing the operator of 



pseudo-gradient with the regular gradient. In the limit l,m function dll^\{—P) can be 
approximated as 



1 



(7rL)V4 



cot - 



9, 



(18) 



where we took into account that ji = 9p. Consequently, the field coefficients a„, for m — ±l can 



be presented in terms of the m - 
that products of the form fl^a„ 



: coefficient ao as a„j = ao[cot9p/2]'". It is clear, therefore, 
in this case are purely real, so that any terms in the force 



proportional to ^m[a*ja,„i] vanish. Since the field strength decays quickly when the particle 
moves out of the equatorial plane 9p = n/2, we only consider = 0^ — 7r/2 <C 1, and expand 

d\l^\{—9p) in terms of 9. However, since for L 3> m, L9 ^ 9 we shall only expand the terms, 
which do not contain L. Keeping linear terms in the expansion of [cot 0^/2]™, and taking into 

account that for L 1 the imaginary part of the Hankel function h^j^\p) is much greater than 
its real part in the region p < L, we obtain the gradient portion of the optical force P'^) = V (u) : 

1 ^ ,. L 



1 



[niikrp)] 



kr. 



-nL{krp)99 



where 



1 

4^ 



^e[a\\EQ\^kVLnL{krp)e 



(19) 



(20) 



The leading term in ac^J^m[{E* ■ V)E] is of the second order in 9 and can be neglected 



and ni^{krp) is the spherical Neumann function 
ad 

compared to the scattering force, F^* 



ac (g) , which takes the form: 
p L 



yl 



1 kr, 



niikrp)^ 



(21) 



where p ~ 2k^a()/3 ^ 1. Thus, in this approximation, the gradient force draws the particle 
toward the resonator radially, and maintains it in the equatorial plane via its 9 component. It 
is important to note that its dependence on the particle's position follows the behavior of the 
Neumann function nL{kr), which is almost exponential in the considered range of parameters. 
At the same time, the azimuthal component of the force is of purely scattering nature and pro- 
portional to the moment density, which can be presented in the simple form (g) — {L/kr)ue,„^, 
where Uem = fio (|2;p) is the electromagnetic energy density (the use of (£ symbol for the field 
emphasizes that it is not affected by the presence of the particle). 

3.3. Effect of the particle-induced modification of the WGM on the optical force 
3.3.1. Modification of the WGM by the particle 

The problem of determination of the electromagnetic field of the coupled resonator-dipole sys- 
tem is analytically tractable and was solved in Ref p6l[27) . The dipole is modeled as a small 
sphere with radius Rp, where kRp ^ 1, and refractive index «p, and the field is found in the form 
of a general expansion in terms VSHs with all l,m and polarizations. The particle is found to 
modify the WGM in two significant ways. First, it creates an additional resonance at frequency 
(Op = (oj^^ + 5cOl with width Fp = + SF^ in addition to the (oj^^ resonance of a single 
sphere, with 5(0l and 5Fi depending only on rp. Second, the steady state of the resonator field 
associated with the cOp resonance is significantly modified compared to the field distribution of 
the initial WGM. Initially isotropic field turns into a highly directional distribution oriented pre- 
dominantly toward the particle. Thus, a displacement of the particle in the ^ direction causes 



the resonator's field to move with the particle (see for details Ref. |26 27 1). In addition, the 
interaction with the particle excites in the resonator WGMs with different I and polarization. 
However, these contributions are small, and can be neglected. In this approximation, the scat- 



tered field of the resonator can again be presented in the form of Eq. ( 14 1, but with expansion 



coefficients, which are no longer given by Eq. ( fTS) . They have the following form 



[yo + i] m=/=±l 

^ L 

Fn 



= - Ff [V^.i(^^p)]'; 5Fl = p\5col\ (23) 



-'e-^<i-0.)>^\ Ffy^_, \ (22) 

where y = {(O — (Op)/Fp. The frequency shift and additional broadening of the resonance 
are Il26]|27): 

^e[a]k ^^(Q)^ 
6neo 

where Vi^iikrp) is the VSH translation coefficient p5) , which arises when the field scattered 
by one sphere is expressed in terms of VSH centered about the other, and is given by 

/3 

VL,n{krp) ^ i{-lf+'^mV2L+lhL{krpy'^ « i-l)'^^/3L/2nLikrp). 
3.3.2. Calculation of the force with particle-modified field 

In this part of the paper we assess the effects of the particle-induced shift of the resonance 
frequency and of the changes in the spatial configuration of the field of the resonator on the 
optical forces exerted by it. To this end we shall analyze the expressions for the force obtained 



by evaluating Eq. ( 10 1 with the field at the location of the particle given by Eq. ( 14 1. The role 
of the pseudo-gradient operator in this equation is to distinguish between field coordinates r 
and particle coordinates Tp even though we calculate the force at the point r = rp. Taking into 
account that dependence on rp is only contained in the expansion coefficients a^, this procedure 



becomes rather trivial: one needs to find the gradient of all respective expressions treating these 
coefficients as constants, and after that equate r = Tp. Calculating the required gradients we 
obtain for the pseudo-gradient component of the force: F'^''-'''' = V(m) 
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The scattering component of the force, F^'', takes in this case the form of 



F(^'-eo|£o| 



3r 



- \h'l'{kr)Y { BJ'm [(fl5(fl_i - fli )] + ^Me [{al{ai + fl_ 1 )] } (27) 



A contribution from the remaining term in Eq. (lOi, which is proportional to ^m(E* • V)E 



remains negligible and will not be considered any further 

One can see that Eq. and Eq. ( [27| ) differ significantly from the respective Eq. ([T9]l and 
Eq. (21 1 obtained under the assumption of the unmodified WGM. Further analysis of the ob- 
tained expression will be performed in two steps. Since it is often assumed that the main effect 
of the particle on the cavity mode consists in changing the resonance frequency, we first separate 
this effect. To achieve this we allow the particle to shift and broaden the resonance according 
to Eq. ( [23| , but will assume that the field coefficients are given by the unmodified Eq. ( [T5] l 
with replacement of (yo + O^^ by (rp/r|^'')(3' + /)^'. In this case the particle can modify the 
amplitude of the resonator's field, but does not change its spatial distribution. Then, one im- 
mediately sees that in Eq. ( [24| , F^''^\ which contains terms proportional to ^m[fl,„fl,„/] with 
m^m' vanishes. For the same reason the 0-component of the scattering force also vanishes. 

In order to present the radial component in the form easily related to expressions used in 
papers of other authors such as those of Ref. pO 13 28|, we introduce the power radiated by 
those modes of the resonator that interact with the particle as. 



and the number of photons in them as 
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The expression for the radial component of the force can be written down now as 

d50}L 



dr„ 



(28) 



which should be contrasted with the expression obtained from the usual derivative of interaction 
energy 

dupoi ^^^ ddcoi dN ^ 



-Nh- 



drp drp 



(29) 



The last terni in this expression is spurious as was first noticed in Ref. p8| , where it was shown 
by direct numerical simulations that the force on a dipole in Fabry-Perot resonator must be 



given by Eq. ( 28 1 rather than by Eq. ( 29 1 



This result is a clear demonstration of the fact that the polarization energy of the particle 
cannot be considered a true potential energy even if one neglects the spatial modification of 
the cavity mode due its interaction with the particle. However, taking this modification into 
account results even in more drastic changes in the optical force yielding a non-zero azimuthal 
component of the pseudo-gradient force and a non-zero polar component of the scattering force. 
The expression for the total force F = F^''*'' +F'^'' in this case is found by using the correct set 



of the field coefficients as defined by Eq. (22i. The radial component of the force does not 



change from Eq. ( |28] l, while its polar and azimuthal components now can be presented as 

r,rf\yl + l){y 



Fe = 2hN — j^^-^ [ScoLil +yyo) + 5rL{yo-y)] (30) 



F^^lhN — ,,,, [5coL{yo-y) + dTL{l+yoy)] (31) 

rprf\yl + l){y^ + l) 

Terms proportional to 5cOl come from the psuedo-gradient force, while STl terms come from 
the scattering force. These expressions demonstrate significant deviation of the force from both 
completely unmodified and frequency-only modified WGM approximations. First, let us note 
that the radial dependence of the force is determined by the factor (y^ + 1 ) ^ ' in addition to 
the Neumann function in 5cOl and STi. The role of this factor can be seen as follows. The 
condition y is satisfied for some = ro at which driving frequency co coincides with the 
particle-induced resonance. If one linearizes y{ri,) about this point as y = (r^ — '"o)/ ('"o), where 

If p is sufficiently small, then in the region where / (ro) can be considered constant, the factor 
ddcoi/drp in iv is also constant. Therefore the spatial profile of the force has Lorentzian shape 
peaked at Vp = ro with width 1 // . Let us recall that in the unmodified WGM approximation 
the magnitude of the force monotonically (essentially exponentially) decreases with rp. 

The azimuthal force is no longer solely due to the scattering contribution. Two differ- 
ent limits are of interest based upon choice of the external driving frequency (0. In the limit 
(0 — >■ (of^ , y5(0i + 5Ti. The magnitude of the pseudo-gradient term exceeds that of the 
scattering term unless y^p- This can only happen for very small values of 5cOi, when the mag- 
nitude of the force is also very small. When y p, the pseudo-gradient contribution exceeds 
the scattering force, and the tangential component can be written as F^j, — {F^^^F^I^'^)/ {ypFp), 



where ^i"' is the scattering force in the unmodified WGM approximation as given by Eq. ill 



If y satisfies 1 <C y < T^l^^ /{pFp), the tangential force exceeds F^^\ When co — (O'p ^ F^" the 
scattering contribution to becomes negligible as well. In this case the force can be written 

~ ^l^^^L^ / P'^p^ which also exceeds the magnitude that would have been obtained in the 
unmodified WGM approximation. These results show that the force propelling the particle in 
the experiments like the one of Ref. p2) is not necessarily of scattering origin and might have 
a pseudo-gradient contribution. The two can be distinguished by their dependence on Oq: while 
the pseudo-gradient force is linear in this parameter, the scattering force is quadratic. 

Even in the range of parameters where the pseudo-gradient contribution to the azimuthal 
component of the force dominates, it remains non-conservative since it imparts net kinetic 



energy to the particle moving along a closed orbit around the resonator. This occurs because 
there is a field gradient which pushes the particle in the (j) direction. When the particle moves to 
a new point, the field re-adjusts so that there is again a field gradient in the direction. Implicit 
in this analysis is the assumption that the particle moves slowly enough to consider the field 
always remaining in the quasi-steady state. Velocity dependent effects can become significant 
when the time scale of particle motion is comparable to the relaxation time of the resonator, 

The relative magnitudes of the force components can be analyzed by comparing {L/rp)5(0L 
to d8(i)i/drp. From the asymptotic expansions for the spherical Nuemann functions and their 



derivatives in the region I, kr < L, we have 11^/11'^^ —cosh(a) \ 29 1 where prime denotes 
differentiation with respect to argument and a is defined by krp = (L+ \ /2)sech{a). Since 
resonances are in the region krp k, L/n, where n is the refractive index of the resonator, we 
have ni/n'^ w — n and thus {d5cOL/drp)/{Ldo}L/rp) « ~2nkrp/L « 2. Thus, assuming we are 
near a particle induced resonance so that < [yj < (9(1), the relative magnitude of the forces 
will be determined by the factor (jq + 1) ' in ^0 and Fq. If the system is driven at a frequency 
close to the ideal Mie frequency, so that yo is of order unity, then the radial and azimuthal forces 
will be of comparable magnitude, while the polar force will be smaller by a factor O.lfyo^ 1 
on the other hand, then the azimuthal force will be smaller than the radial force by a factor 
1 /yo, while the polar force is smaller by a factor Oy/yo- 

It can also be seen that the scattering contribution to the azimuthal and polar forces is in 
general smaller than the pseudo-gradient contribution due to the fact that STl/ScOl = p ^1 
(where p — 2A;^ao/3). In the limit where both the driving frequency and particle induced res- 
onance frequencies are very close to the ideal Mie resonance, so that y,yo — > 0, the scattering 
contribution to becomes appreciable, while it vanishes in Fq. This is to be expected given 
that y,yQ -^Q is the limit where the particle induced modification of the cavity mode becomes 
vanishingly small, and accordingly the behavior of the forces approaches that of their unmodi- 



fied forms of Eq. ( 19 1 and Eq. (21 1. This is likely the regime encountered in the experiments of 
Ref. 

The results of the calculation of the force within the pseudo-gradient approach can be com- 
pared with calculations carried out by integrating the Maxwell stress tensor over a surface of 
the particle, which is assumed to have a small, but finite size. For a field represented by a VSH 
expansion, the stress tensor integral over a spherical region can be performed analytically and 
the force given in terms of the VSH expansion coefficients |30|. For the present case these co- 
efficients are given in Ref. p7) while full details of the calculations of the force can be found 



in Ref. |21 1. In the large L limit, the forces obtained agree exactly with those calculated from 
the pseudo-gradient approach validating the latter. 

It is interesting to note that the large L limit of the stress tensor calculations is necessary to 
maintain consistency with the assumed point-like nature of the particle in the pseudo-gradient 
approach. To see this, note that for a given resonator of radius R and refractive index n, the 
lowest order approximation to the resonant frequency is the geometric optical condition nkR w 
L. At the same time, a point dipole is defined by the the limit Rp^O with electromagnetic size 
parameter p = kRp kept constant. Combining these two conditions we have nRp = LRp, which 
implies that taking — > requires that L — > 00. 

4. Conclusion 

We have presented here a generalization of the theory of optical forces on a dipole for the 
case when it interacts with the electromagnetic field of an optical cavity. The traditional gra- 
dient/scattering paradigm is shown to be invalid when the dipole can modify the source of the 
field. In particular, all vector components of the force are found to be non-conservative and. 



consequently, no component can be derived from a gradient of electromagnetic polarization 
energy. We have further shown that, when the particle-induced modification of the resonator 
field is taken into account, the force in the direction of the energy and momentum flux of the 
wave cannot be interpreted as a 'scattering' force. In place of the gradient/scattering paradigm 
we have proposed a pseudo-gradient framework which is conceptually simpler and computa- 
tionally more efficient than the exact Maxwell stress tensor approach. Using an example of 
a small dielectric particle interacting with whispering-gallery-modes of a spherical resonator 
we demonstrated that the suggested pseudo-gradient formalism reproduces all results of the 
calculations based on the Maxwell stress tensor. 

The results of this work have important implications for the quantum theory of optomechani- 
cal interaction, which is commonly based on the assumed potential nature of the gradient force. 
These results are also of importance for proposed optofluidic sensors which rely on a tangential 
force to drive the particle in orbit around the resonator. 



